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Abstract 

We extend a recent sum rule calculation for inelastic quarkonium-hadron inter- 
actions to realistic parton distribution functions; we also include finite target-mass 
corrections. Both modifications are shown to have no significant effect on the re- 
sulting cross section behaviour but the performed analysis gives useful insights on 
the sum-rule approach in general. 
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In a recent study the cross section <7jv$(s) for inelastic quarkonium-nucleon collisions 
was calculated from sum rules established on the basis of the operator-product expansion 
0-0; here $ denotes a quarkonium ground state (J/0 or T) and s the squared center of 
mass collision energy. In order to simplify the calculations a number of approximations 
was made in this investigation. The aim of the present work is to obtain 07^(5) avoiding 
such approximations and to clarify some aspects of the sum rule approach in general. 

In the formalism developed to describe the interaction of hadrons with heavy qq mesons 
@) 0> HI 7 the large quark mass vhq allows one to write down the corresponding Born 
amplitude in the form of operator product expansion (OPE). The effects governed by 
short-distance physics, including the structure of heavy quarkonium state, factorize into 
coefficient functions, C n , while the thorough dependence on soft structure of hadron is 
determined by low-energy matrix elements of gauge-field local operators, O n , renormalized 
at quarkonium mass scale M$, 
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■M-Born — 7 y C n {O n ). 



n=2,4,. 



The operators entering (1) are ordered by their dimensions. If the characteristic distances 
of internal qq dynamics are small enough as compared to scale of QCD-vacuum fluctua- 
tions inside the hadron, only the terms of lowest dimension contribute. In the framework 
of the multipole expansion formalism J7J, [[5], the lowest dimension terms correspond to 
the dipole field contribution. For the sake of simplicity, we have dropped in (1) spin 
indices which otherwise mark the amplitudes describing the forward scattering of heavy 
quarkonia with different polarizations. If we restrict ourselves to the case of spin-averaged 
interactions, then Lorentz indices may appear only in the following combination: 

® n = ~]\/[n~^Vl ■ ■ ■ Knn@ni...n n l (2) 



where 



qG = -n-2 



G aiU V^ 2 . . . ^V-i^Wj (symmetrized) - traces. (3) 

Here is the quarkonium momentum, G is the gluon field operator and T> means the 
covariant derivative. 

The target matrix elements of these operators are completely specified by pheno- 
menological parameters A n of dimension defined as twist and the spin-averaged tensor 
structure, 

(P\°Z-^ \P) = A n G° W • • • P 11 " ~ traces ) . ( 4 ) 

thus forming into traceless, symmetrical rank-n tensors. Substituting (4) and (2) into (1) 
we come to schematic expression 

m, p = £ CA ^^ n^*, (5) 

n=2,4,... 1V1 $> 

where for convenience we used the notation IP 1 "^" = (p^ 1 . . .p^ n — traces). 

At this point, we comment briefly on the applicability of the analysis to actual physical 
states. For sufficiently heavy constituent quarks, the state $ is localized at characteristic 
distances r which are sufficiently small compared to non-perturbative hadronic scale Aq^ d 
to treat the system perturbatively. On the other hand, it satisfies r 3> mn , so that one can 
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describe the internal dynamics nonrelativistically. Although the estimates show that 
these requirements are reliably satisfied only for quark masses exceeding 25 GeV, one finds 
that a Coulomb-like potential approach provides a satisfactory description of quarkonium 
spectroscopy (see, e.g. 0). Hence we may also expect a one-gluon exchange picture to be 
relevant for J/ip and T internal dynamics. We thus consider $ as SU(3)-Coulombic bound 
state characterized by a "Bohr" radius ao and the corresponding "Rydberg" energy e : 

ao = o > e o = -H m Qi (6) 
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where a s is the gauge coupling constant, evaluated at a scale e . Then the direct calcu- 
lations of |H in SU(N) gauge theory give 

for the Wilson coefficients C n of lS'-quarkonium state in the leading order in 1/N 2 (anal- 
ogous expressions for the P-states can be found in Ref. [[10]]). 

We now insert the coefficients (7) into Eq. (5). If we neglect the trace terms in Eq. 
(4), or, equivalently, the corrections of order of m^/eg, we recover the sum rules used in 

as the basis for calculating <tn$(s), 



dy y n 2 \Jl- y 2 a N ^(m N /y) = I(n) j dx x n 2 g(x,Q 2 
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I(n)=2^ 2 ^ -±-41 -1- -L.. (9) 



/o " J 
with I in) given by 

16\ 2 r (n + I) /j^x 

3 / T(n + 5) V3aJ 

One obvious problem in applying these sum rules to J/ip interactions is that ratio m 2 N /e^ ~ 
2.1 is not actually small and hence cannot be ignored. In other words, the trace terms 
entering the definition of twist-two operators (4) in the Wilson's ordering scheme must 
be included explicitly producing the corresponding changes in dispersion sum rules (8). 
The tracelessness of the tensor IP 1 '^™ () means that 

//^Ji'' --'' --'' = 0. (10) 

The most general structure of IP 1 '"^" is well known 

n/2 / -\ j j n-2j 

n w ~"» = E(-!) ii! ^#^ E T^~ui^~i>- (ii) 

j=0 " " all permut. 

where the second sum runs over the n\/(n — 2j)!(2j)! terms of all possible permutations. 
Introducing A = (pK/M§) we can rewrite (5) as 

av = tj^lyiy^^^ (fj gp . (12) 
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Changing the summation index n — > n — 2j and making simple rearrangements of the 
coefficients, one obtains 
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n=0,2,... \ e °/ j=0 



[n + j)\ ( m 2 N V 



(13) 



We now take the explicit expressions (7) for d n (the term proportional to do corresponds 
to disconnected diagrams and does not enter the amplitude of quarkonium interactions in 



the OPE formulation) and use the definition of A n as Mellin transforms [12] of the gluon 
distribution in a proton, evaluated at the scale Q 2 = 

A n = dx x n ~ 2 g(x,Q 2 = el). 
Jo 

This leads to the final expression 



Mfr = 40rao e o (y) 



f dx £ x n - 2 (±) n g( 

J ° n=2,4,.. \ e 0/ 



x,e )x 



r(n + f) 

T(n + 5) 



m 



5 n 7 n (5 + n) 
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Introducing the variable y = m N /\, we get 



(14) 



(15) 



dyy^il-y^aM^N/y) 



Kn) 



n 



dxx n 2 g(x, el) 3 F 2 [^ + ^,- + ^,1 + 71; 
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n 



(5 + n) n_ 
2 ' 2' 



(16) 



Comparing Eqs. (8) and (16), we see that the inclusion of finite mass corrections effectively 
modifies the gluon distribution function - in the r.h.s. one should integrate including 
confluent hypergeometric function 3-F 2 (. . .) as a weight. In Eq. (16), n = 2,3, so that 
in order to determine the cross section, it appears that we have to solve an infinite set 
of equations. However, usually the application of the sum rules to deep-inelastic lepton- 
hadron scattering [TB| and heavy meson photoproduction |5] had been limited to the 
values of n < 5 — 6 [1-5]: it was expected that perturbative calculation of the coefficient 
functions breaks down for large n. The use of Eq.(16) can effectively extend the range 
of applicability of the sum rule approach, and the problem of relevant n-interval should 
be re-examined. Indeed, the additional weight function 3 F 2 in the r.h.s. of (16) decreases 
for x — >• 1 faster for larger n, making the sum rules more sensitive to the behaviour of 
parton distribution function (PDF) at lower values of x. This makes the reliability region 
for gluonic distribution function moments effectively wider. In what follows we study the 
sum rules for values n < 10. 

It is evident from Eq. (16) that the behaviour of the cross section at high energy is 
particularly sensitive to the small x region, while the threshold behaviour probes the large 
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x region. New data from deep inelastic scattering have led to different PDF parametriza- 
tions; hence one has to check how well the sum rules, including target mass corrections 
as well, are satisfied for different PDF choices. In the simple schematic form 

g(x) = const. (A; + 1)(1 -a;) 4 (17) 

was used. Here, we have considered in addition the MRS DJ and the new MRS H [ |I4| 
parametrizations; the latter takes into account the small- x behavior observed at HERA 
TE \. In all cases, the left hand and right hand sides of Eq. (16) agree for 2 < n < 10 
within better than 1%. The sum rules are thus well satisfied for all three PDF forms used. 

Next we want to check what effect the finite mass corrections have. In Fig. 1, we 
therefore show the ratio of the cross section including the target mass corrections (i.e. 
the solution of Eq.(16)) to the same calculation without them (from Eq.(8)). It is seen 
that the inclusion of target mass corrections changes the resulting cross section above the 
threshold (y/s > 5 GeV) by less than a factor three. 

Finally we want to see if and how the cross sections obtained as solution of Eq. (16) 
vary for different PDF forms; we have therefore calculated <7jv$(s) also for MRS D'_ and 
MRS H. The resulting cross sections, divided by the corresponding form obtained using 
PDF (17), are shown in Fig. 2; here all cross sections include target mass corrections. 
The increase of the two MRS forms at high energy, relative to that using Eq. (17), is due 
to the small x increase in the more realistic PDF's. 

In summary, we would like to stress that including finite target mass corrections jus- 
tifies the use of OPE close to the threshold and meanwhile does not change significantly 
the behaviour of quarkonium-hadron cross section at high energies. Besides, target mass 
corrections are important in providing the basis for future investigations of the role of 
higher terms of the multipole expansion and for a description of the entire x region by 
the sum-rule method. 
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Figure Captions: 

Fig. 1: The ratio of <7j/$_jv(s), including target mass corrections, to the corresponding 
form without these corrections; both forms are calculated with PDF (17). 

Fig. 2: The ratio of aj/^ N (s), from MRS D-' PDF (solid line) and from MRS H (dashed 
line), to the cross section obtained from PDF (17); all cross sections include target mass 
corrections. 
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